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ABSTRACT 

We compute new chemical profiles for the core and envelope of white dwarfs appropriate for pulsa- 
tional studies of ZZ Ccti stars. These profiles are extracted from the complete evolution of progenitor 
stars, evolved through the main sequence and the thermally-pulsing asymptotic giant branch (AGB) 
stages, and from time-dependent element diffusion during white dwarf evolution. We discuss the 
importance of the initial-final mass relationship for the white dwarf carbon-oxygen composition. In 
particular, we find that the central oxygen abundance may be underestimated by about 15% if the 
white dwarf mass is assumed to be the hydrogen-free core mass before the first thermal pulse. We 
also discuss the importance for the chemical profiles expected in the outermost layers of ZZ Ceti stars 
of the computation of the thermally-pulsing AGB phase and of the phase in which element diffusion 
is relevant. We find a strong dependence of the outer layer chemical stratification on the stellar mass. 
In particular, in the less massive models, the double-layered structure in the helium layer built up 
during the thermally-pulsing AGB phase is not removed by diffusion by the time the ZZ Ceti stage 
is reached. Finally, we perform adiabatic pulsation calculations and discuss the implications of our 
new chemical profiles for the pulsational properties of ZZ Ceti stars. We find that the whole g— mode 
period spectrum and the mode-trapping properties of these pulsating white dwarfs as derived from 
our new chemical profiles are substantially different from those based on chemical profiles widely used 
in existing asteroseismological studies. Thus, we expect the asteroseismological models derived from 
our chemical profiles to be significantly different from those found thus far. 

Subject headings: stars — pulsations — stars: interiors — stars: evolution — stars: white dwarfs 



I. INTRODUCTION 

Pulsating DA (H-rich atmospheres) white dwarfs, also 
called ZZ Ceti or DAV stars, are the most numerous class 
of degenerate pulsators, with over 143 members known 
today (Winget & Kepler 2008). Since the discovery of 
the first ZZ Ceti star, HL Tau 76, by Landoff (1968), 
there has been a continuous effort to model the interior 
of these variable stars. ZZ Ceti stars are found within 
a very narrow strip of effective temperatures (10 500 K 
^ ^eff ^ 12 500 K). They are characterized by multiperi- 
odic brightness variations of up to 0.30 mag caused by 
spheroidal, non-radial g-modes of low degree {I < 2) with 
periods between 100 and 1200 s. The driving mechanism 
thought to excite the pulsations near the blue edge of 
the instability strip is the k — 7 mechanism that takes 
place in the hydrogen partial ionization zone (Dolez & 
Vauclair 1981; Dziembowski & Koester 1981; Winget et 
al. 1982). Also, the "convective driving" mechanism has 
been proposed — first by Brickhill (1991) and later re- 
examined by Goldreich & Wu (1999). It appears to be 
the responsible of mode driving once a thick convection 
zone has developed at the stellar surface. 

The comparison of the observed pulsation periods in 
white dwarfs and the periods computed for appropri- 
ate theoretical models (white dwarf asteroseismology) al- 
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lows to infer details of their origin, internal structure and 
evolution (Winget & Kepler 2008; Fontaine & Brassard 
2008). In particular, the stellar mass, the thickness of the 
outer envelopes, the core chemical composition, magnetic 
fields and rotation rates can be determined from the ob- 
served periods. In addition, the asteroseismology of ZZ 
Ceti stars is a valuable tool for studying axions (Isern et 
al. 1992; Corsico et al. 2001; Bischoff-Kim et al. 2008; 
Isern et al. 2010) and crystallization (Montgomery et 
al. 1999; Corsico et al. 2004, 2005; Metcalfe et al. 2004; 
Kanaanetal. 2005). Finally, the temporal changes in the 
observed periods can help detect planets orbiting around 
white dwarfs (Mullally et al. 2008). 

The first published complete set of DA white dwarf 
models suitable for asteroseismology was that of Tassoul 
et al. (1990). A large parameter space was explored in 
such a monumental study, and for a long time (since the 
early eighties) this set of models represented the state- 
of-the-art in the area. The pulsation properties of these 
models were thoroughly explored in a series of important 
papers by Brassard et al. (1991, 1992a, 1992b). As im- 
portant as these models were at that time, they suffer 
from a number of shortcomings. For instance, the core 
of the models is made of pure carbon, while stellar evo- 
lution calculations indicate that cores of typical white 
dwarfs are made of a mixture of carbon and oxygen. 
Also, the carbon/helium (C/Hc) and helium/hydrogen 



(He/H) chemical interfaces are modeled on the basis of 
the assumption of the diffusive equilibrium in the "trace 
element approximation" , an approach that involves a 
quasi-discontinuity in the chemical profiles at the tran- 
sition regions which, in turn, leads to peaked features 
in the Brunt- Vaisala frequency and exaggerated mode- 
trapping effects (Corsico et al. 2002a, 2002b). These 
models were employed for asteroseismological inferences 
of the DAVs G 226-29 (Fontaine et al. 1992) and GD 
154 (Pfeiffer et al. 1996). More recently, Pech et al. 
(2006) and Pech & Vauclair (2006) have presented as- 
teroseismological analysis on HL Tau 76 and G 185—32, 
respectively, by employing similar DA white dwarf mod- 
els, although with updated input physics. 

The models of Bradley (1996) constituted a substan- 
tial improvement in the field. These models have carbon- 
oxygen cores in varying proportions, and the C/He and 
He/H chemical interfaces are more realistic. Perhaps the 
most severe shortcoming of these models is the (unrealis- 
tic) ramp-like shape of the core carbon-oxygen chemical 
profiles. These DA models were the basis of the very im- 
portant asteroseismological studies on the DAVs G 29- 
38 (Bradley & Kleinman 1997), G 117-B15A and R 548 
(Bradley 1998), GD 165 and L 19-2 Bradley (2001), and 
G 185-32 (Bradley 2006). 

The next step in improving the modeling of DAVs 
was given by Corsico et al. (2002b) and Benvenuto et 
al. (2002a), who employed evolutionary models char- 
acterized by He/H chemical interfaces resulting from a 
time-dependent element diffusion treatment (Althaus & 
Benvenuto 2000), and the carbon-oxygen core chemical 
structure extracted from the evolutionary computations 
of Salaris et al. (1997). The use of very smooth outer 
chemical interfaces, as shaped by chemical diffusion, re- 
vealed that the use of the trace element approximation 
turns out to be inappropriate to model the shape of the 
chemical interfaces in a DA white dwarf. This grid of 
models was employed in an asteroseismological study of 
G 117-B15A (Benvenuto et al. 2002a). For these se- 
quences, the starting configurations for the white dwarf 
evolution were obtained through an artificial procedure, 
and not as result of evolutionary computations of the 
progenitor stars. 

Recently, Castanheira & Kepler (2008, 2009) have car- 
ried out an extensive asteroseismological study of DAVs 
by employing DA white dwarf models similar to those 
of Bradley (1996), but with a simplified treatment of 
the core chemical structure, by somewhat arbitrarily 
fixing the central abundances to 50 % oxygen and 50 
% carbon. The He/H chemical interfaces adopted for 
these models are a parametrization of the realistic chem- 
ical profiles resulting from time-dependent element dif- 
fusion (Althaus et al. 2003). The study includes the 
"classical" DAVs and also the recently discovered SDSS 
DAVs. In total, 83 ZZ Ceti stars are analyzed. An im- 
portant result of these studies is that the thickness of 
the H envelopes inferred from asteroseismology is in the 
range 10^"* > Afn/Af* J^ 10^^", with a mean value of 
Afn/Af* = 5 X 10^^. This suggests that an important 
fraction of DAs characterized by envelopes substantially 
thinner than predicted by the standard evolution theory 
could exist, with the consequent important implications 
for the theory of white dwarf formation. However, these 



results are preliminary and do not include the possible 
effects of realistic carbon-oxygen profiles on the astero- 
seismological fits. 

Almost simultaneously with the study of Castanheira 
& Kepler (2008, 2009),'Bischoff-Kim et al. (2008) per- 
formed a new asteroseismological study on G 117— B15A 
and R 548 by employing DA white dwarf models simi- 
lar to those employed by Castanheira & Kepler (2008, 
2009), but incorporating realistic core chemical profiles 
according to Salaris et al. (1997). The results of this 
work are in reasonable agreement with previous studies 
on these ZZ Ceti stars. However, the mass and effec- 
tive temperatures found by Bischoff-Kim et al. (2008) 
for G 117— B15A are rather high (especially the mass, at 
0.66 M0). Recently, Bognar et al. (2009) have employed 
the same asteroseismological modeling to study the pul- 
sations of the ZZ Ceti star KUV 02464-1-3239. Finally, 
Bischoff-Kim (2009) presented the results of an astero- 
seismological analysis of two DAVs with rich pulsation 
spectrum, G 38—29 and R 808 based on similar mod- 
els, with parametrized, smooth ramp-like core profiles. 
These models are able to reproduce the observed period 
spectra reasonably well, though some assumptions about 
the 771 and £ identification of modes were made. 

White dwarf stellar models with realistic chemical pro- 
files are crucial to correctly assess the adiabatic period 
spectrum and mode-trapping properties of the DAVs, 
which lies at the core of white dwarf asteroseismology 
(Brassard et al. 1992a; Bradley 1996; Corsico et al. 
2002a). In this paper, we compute for the first time 
consistent chemical profiles for both the core and the 
envelope of white dwarfs with various stellar masses ap- 
propriate for detailed asteroseismological fits of ZZ Ceti 
stars. These chemical profiles are extracted from the 
full and complete evolution of progenitor stars from the 
zero age main sequence, to the thermally-pulsing and 
mass-loss phases on the asymptotic giant branch (AGB), 
and from time-dependent element diffusion predictions 
during the white dwarf stagcQ- These profiles will be 
valuable in conducting future asteroseismological studies 
of ZZ Ceti stars that intend to include realistic chemi- 
cal profiles throughout the interior of white dwarfs. To 
assess the impact of these new profiles on the theoreti- 
cal pulsational inferences, we perfom adiabatic pulsation 
computations and compare the resulting periods with 
the pulsational inferences based on the most widely used 
chemical profiles in existing asteroseismological fits. 

The paper is organized as follows. In Sect. [21 we pro- 
vide a description of the input physics assumed in the 
evolutionary calculations of relevance for the chemical 
composition. In Sect. |3] we discuss the importance of 
the initial-final mass relationship for the expected white 
dwarf carbon-oxygen composition. The resulting chemi- 
cal profiles are described at some length in Sect. 21 The 
implications of our new chemical profiles for the pulsa- 
tional properties of ZZ Ceti stars are discussed in Sect. 
[5l We conclude in Sect. [Hby summarizing our findings. 

2. INPUT PHYSICS 



^ The chemical profiles for the core and envelope of 
our m odels appropriate for ZZ Ceti stars are available at 
|http : //www . f caglp . unlp . edu . ar/evolgroup 



TABLE 1 

Initial and final stellar mass (in solar units), and the 

CENTRAL OXYGEN ABUNDANCE (mASS FRACTION) LEFT AFTER 

CORE HELIUM BURNING, AND AFTER RAYLEIGH- TAYLOR 

REHOMOGENIZATION. THE PROGENITOR. METALLICITY IS Z = 0.01. 



JW'ZAMS 


MwD 


Xo(CHB) 


Xo(RT) 


1.00 


0.5249 


0.702 


0.788 


1.50 


0.5701 


0.680 


0.686 


1.75 


0.5932 


0.699 


0.704 


2.00 


0.6096 


0.716 


0.723 


2.25 


0.6323 


0.747 


0.755 


2.50 


0.6598 


0.722 


0.730 


3.00 


0.7051 


0.658 


0.661 


3.50 


0.7670 


0.649 


0.655 


4.00 


0.8373 


0.635 


0.641 


5.00 


0.8779 


0.615 


0.620 



The chemical profiles presented in this paper have been 
extracted from full and complete evolutionary calcula- 
tions for both the white dwarfs and their progenitor stars, 
by using an updated version of the LPCODE stellar evo- 
lutionary code — see Althaus et al. (2005) and refer- 
ences therein. This code has recently been employed 
to study different aspects of the evolution of low-mass 
stars, such as the formation and evolution of DA white 
dwarfs (Rcnedo et al. 2010), H-deficicnt white dwarfs, 
PG 1159 and extreme horizontal branch stars (Althaus 
et al. 2005; Miller Bertolami & Althaus 2006; Miller 
Bertolami et al. 2008; Althaus et al. 2009a), as well as 
the evolution of He-core white dwarfs with high metallic- 
ity progenitors (Althaus et al. 2009b). It has also been 
used to study the initial-final-mass relation in Salaris et 
al. (2009), where a test and comparison of LPCODE with 
other evolutionary codes has also been made. Details of 
LPCODE can be found in these works. In what follows, we 
comment on the main input physics that are relevant for 
his work. We assume the metallicity of progenitor stars 
to be Z = 0.01. 

The LPCODE evolutionary code considers a simultane- 
ous treatment of non-instantaneous mixing (and extra- 
mixing if present) and burning of elements (Althaus et 
al. 2003). The nuclear network accounts explicitly for 
the following elements: ^H, ^H, ^He, ^He, ^Li, ^Be,' ^^C, 
13c, "N, i^N, 160, "O, 1^0, i^F, 20Ne and 22Ne, to- 
gether with 34 thermonuclear reaction rates of the pp- 
chains, CNO bi-cycle, helium burning, and carbon ig- 
nition that are identical to those described in Althaus 



et al. (2005), with the exception of i^C -|- p^- 



13 



N 



13f 



7 -^ ""^C + e"*" -I- j/q and ^C(j>,j) N, which were taken 
from Angulo et al. (1999). The 120(0,7)1^0 reaction 
rate, of special relevance for the carbon-oxygen stratifi- 
cation of the emerging white dwarf, was taken from An- 
gulo et al. (1999). The standard mixing length theory 
for convection — with the free parameter a = 1.61 — 
has been adopted. With this value, the present luminos- 
ity and effective temperature of the Sun, at an age of 
4 570 Myr, are reproduced by LPCODE when Z = 0.0164 
and X = 0.714 are adopted — in agreement with the 
Z/X value of Grevesse & Sauval (1998). 

Except for the evolutionary stages corresponding to the 
thermally-pulsing AGB phase, we have considered the 
occurrence of extra-mixing episodes beyond each convec- 
tive boundary following the prescription of Herwig et al. 
(1997). We treated extra-mixing as a diffusion process by 



assuming that mixing velocities decay exponentially be- 
yond each convective boundary. Specifically, we assume a 
diffusion coefficient given by Dqv = Dq exp{—2z/fHp), 
where Hp is the pressure scale height at the convective 
boundary, Dq is the diffusion coefficient of unstable re- 
gions close to the convective boundary, and z is the ge- 
ometric distance from the edge of the convective bound- 
ary (Herwig et al. 1997). We adopted / = 0.016 in 
all our sequences, a value inferred from the width of the 
upper main sequence. The occurrence of extra-mixing 
episodes during core helium burning largely determines 
the final chemical composition of the white dwarf core. In 
this sense, our treatment of time-dependent extra-mixing 
episodes predicts a core chemical stratification similar to 
that predicted by the phenomenon of "semiconvection" 
during central helium burning, which naturally yields the 
growth of the convective core (Straniero et al. 2003). Fi- 
nally, the breathing pulse instability occurring towards 
the end of core helium burning was suppressed — see 
Straniero et al. (2003) for a discussion of this point. In 
our simulations, breathing pulses have been suppressed 
by gradually decreasing the parameter / from the mo- 
ment the helium convective core starts to receed (which 
occurs once the helium abundance in the core decreases 
below « 0.05 — 0.1). At this stage, the gradual suppres- 
sion of extra-mixing toward the end of core helium burn- 
ing bears no consequences for the final chemical stratifi- 
cation. 

Other physical ingredients considered in LPCODE are 
the radiative opacities from the OPAL project (Iglesias & 
Rogers 1996), including carbon- and oxygen-rich compo- 
sition, supplemented at low temperatures with the molec- 
ular opacities of Alexander & Ferguson (1994). During 
the white dwarf regime, the metal mass fraction Z in 
the envelope is not assumed to be fixed. Instead, it is 
specified consistently acording to the prediction of ele- 
ment diffusion. To account for this, we have considered 
radiative opacities tables from OPAL for arbitrary metal- 
licities. For effective temperatures less than 10,000 K we 
have included the effects of molecular opacitiy by assum- 
ing pure hydrogen composition from the computations 
of Marigo & Aringer (2009). This assumption is jus- 
tified because element diffusion leads to pure hydrogen 
envelopes in cool white dwarfs. The conductive opacities 
are those of Cassisi et al. (2007), and the neutrino emis- 
sion rates are taken from Itoh et al. (1996) and Haft et 
al. (1994). For the high density regime characteristics of 
white dwarfs, we have used the equation of state of Scg- 
retain et al. (1994), which accounts for all the important 
contributions for both the liquid and solid phases — see 
Althaus et al. (2007) and references therein. 

In this study, we have considered the distinct physical 
processes that are responsible for changes in the chemical 
abundance distribution during white dwarf evolution. In 
particular, element diffusion strongly modifies the chemi- 
cal composition profile throughout their outer layers. As 
a result of diffusion processes, our sequences developed 
pure hydrogen envelopes, the thickness of which gradu- 
ally increases as evolution proceeds. We have considered 
gravitational settling as well as thermal and chemical 
diffusion — but not radiative levitation, which is rele- 
vant only for the hottest and brightest post-AGB and 
early white-dwarf cooling stages for determining the sur- 
face composition — of iH, ^He, ^He, i^C, "C, I'lN and 
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Fig. 1. — Central oxygen abundance (mass fraction) left after 
core helium burning in terms of both the hydrogen- free core mass 
right before the first thermal pulse (upper panel) and the initial 
stellar mass on the ZAMS (lower panel). The solid lines show our 
results while the dashed lines show those of Salaris et al. (1997) 

^^O, see Althaus et al. (2003) for details. Our treat- 
ment of time-dependent diffusion is based on the multi- 
component gas treatment presented in Burgers (1969). 
In LPCODE, diffusion becomes operative once the wind 
limit is reached at high effective temperatures (Unglaub 
& Bues 2000). In addition, abundance changes result- 
ing from residual nuclear burning — mostly during the 
hot stages of white dwarf evolution — have been taken 
into account in our simulations. Finally, we considered 
the chemical rehomogenization of the inner carbon- oxy- 
gen profile induced by Rayleigh- Taylor (RT) instabilities 
following Salaris et al. (1997). These instabilities arise 
because of the positive molecular weight gradients that 
remain above the flat chemical profile left by convection 
during core helium burning. 

3. THE IMPORTANCE OF THE INITIAL-FINAL MASS 
RELATIONSHIP 

As mentioned, chemical profiles appropriate for DA 
white dwarfs have been derived from the full evolutionary 
calculations of progenitor stars for solar metallicity. To 
this end, the complete evolution of ten evolutionary se- 
quences with initial stellar mass ranging from 1 to 5 Mq 
were computed from the ZAMS through the thermally- 
pulsing and mass-loss phases on the AGB and finally to 
the domain of planetary nebulae. In Table [1] we list the 
stellar masses of the resulting white dwarfs, together with 
the inital masses of the progenitor stars on the ZAMS. 
Also listed in Table [1] is the central oxygen abundance 
both at the end of core He burning and after chemical 
rehomogenization by Rayleigh- Taylor instabilities. 

We mention that extra-mixing episodes were disre- 
garded during the thermally-pulsing AGB phase. In par- 
ticular, a strong reduction (a value of / much smaller 
than 0.016) of extra-mixing episodes at the base of the 
pulse-driven convection zone is supported by simula- 
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Fig. 2. — Initial-final mass relationships: The final mass given by 
the hydrogen-free core mass is depicted in terms of the initial mass 
of the progenitor star. In addition to the observational data from 
open clusters (Catalan et al. 2008), upwards triangles, we show 
the theoretical predictions given by our calculations at the end of 
the TP-AGB phase and before the first thermal pulse (circles and 
squares, respectively). Also shown are the predictions of Weiss & 
Ferguson (2009) and Salaris et al. (1997) before the first thermal 
pulse (stars and downwards triangles, respectively). 

tions of the s— process abundance patterns (Lugaro et 
al. 2003) and, more recently, by observational inferences 
of the initial- final mass relation (Salaris et al. 2009). As 
a result, it is expected that the mass of the hydrogen- 
free core of our sequences gradually grows as evolution 
proceeds through the thermally-pulsing AGB. This is be- 
cause a strong reduction or suppression of extra-mixing 
at the base of the pulse-driven convection zone strongly 
inhibits the occurrence of third dredge-up, thus favor- 
ing the growth of the hydrogen-free core. We considered 
mass-loss episodes during the core helium burning stage 
and on the red giant branch following Schroder & Cuntz 
(2005). During the AGB and thermally-pulsing AGB 
phases, we adopted the maximum mass loss rate between 
the prescription of Schroder & Cuntz (2005) and that of 
Vassiliadis & Wood (1993). In the case of a strong re- 
duction of third dredge-up, as occurred in our sequences, 
mass loss plays a major role in determining the final mass 
of the hydrogen-free core at the end of the TP-AGB evo- 
lution, and thus the initial-final mass relation (Weiss & 
Ferguson 2009). 

We begin by examining Fig. [T] which displays the cen- 
tral oxygen abundance by mass fraction left after core 
helium burning. The upper panel shows the predicted 
central oxygen abundance in terms of the hydrogen-free 
core mass right before the first thermal pulse, while the 
lower panel shows this quantity in terms of the initial 
stellar mass on the ZAMS. The predictions of our cal- 
culations — solid lines — are compared with those of 
Salaris et al. (1997) — dashed lines. Note the qual- 
itatively good agreement between both sets of calcula- 
tions. We recall that the final carbon-oxygen stratifica- 
tion of the emerging white dwarf depends on both the 
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Fig. 3. — Inner carbon and oxygen abundance by mass for the 
0.63 Mq white dwarf resulting from two different progenitors that 
lead to the same white dwarf mass. The chemical profiles corre- 
spond to progenitor stars with initial stellar mass of 3 and 2.25 Mq 
(thick and thin lines, respectively), and at stages before chemical 
rehomogenization by Rayleigh- Taylor instability. 

efficiency of the ^^C(a,7)^^0 reaction rate and tfie oc- 
currence of extra-mixing episodes toward tfie late stage 
of core helium burning. In particular, the systematically 
lower oxygen abundances of our models are due mostly to 
our use of the cross section for the ^^C(a, 7)^^0 reaction 
rate given by Angulo et al. (1997), which is smaller than 
the rate of Caughlan et al (1985) adopted by Salaris et 
al. (1997). Note that both sets of calculations predict a 
maximum in the central oxygen abundance at an initial 
mass of M « 2.5Mq. 

A careful computation of the evolutionary stages dur- 
ing the thermally-pulsing AGB and the resulting initial- 
final mass relationship is an important aspect concerning 
the final carbon-oxygen composition of the white dwarf 
core. This can be seen by inspecting Fig. [5J where vari- 
ous theoretical initial-final mass relationships, giving the 
mass of the hydrogen-free core in terms of the ZAMS 
mass of the progenitor, are plotted. The results shown 
in this figure include the predictions of our full evolution- 
ary calculations at the end of the thermally-pulsing AGB 
phase and at the beginning of the first thermal pulse 
(circles and squares, respectively). Our relationships are 
compared with those of Salaris et al. (1997) and Weiss 
& Fergusson (2009), both at the beginning of the first 
thermal pulse. We also show the semi-empirical initial- 
final mass relationship of Catalan et al. (2008), based on 
white dwarfs in open clusters and in common proper mo- 
tions pairs with metallicities close to Z = 0.01, the metal- 
licity we assume for the progenitor stars of our sequences. 
In view of the discussion of the preceeding paragraph, 
note the increase in the mass of the hydrogen-free core 
during the thermally-pulsing AGB stage. As a result, the 
initial-final mass relationship by the end of the thermally- 
pulsing AGB becomes markedly different from that de- 
termined by the mass of the hydrogen-free core before 
the first thermal pulse. For the carbon-oxygen composi- 
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Fig. 4. — Inner carbon and oxygen abundance by mass for the 
0.525, 0.609, and 0.878 Mq white dwarf models at Tefi ^ 12, 000/C, 
and after chemical rehomogenization by Rayleigh- Taylor instabili- 
ties. 

tion expected in a white dwarf, this is an important issue. 
Indeed, if the mass of the hydrogen- free core is assumed 
to be essentially the mass of the resulting white dwarf, it 
is clear that a white dwarf with a given mass may corre- 
spond to different progenitor stellar masses depending on 
the assumed initial-final mass relationship. For instance, 
from our theoretical initial-to-final mass relationships, a 
white dwarf with 0.63 Mq would correspond to a progen- 
itor star with a stellar mass of 2.25 Mq if the initial- final 
mass relationship is assessed at an advanced stage in the 
thermally-pulsing AGB phase, or 3 Mq if the white dwarf 
mass is assumed to be the hydrogen-deficient core mass 
right before the first thermal pulse. In particular, we 
note that Salaris et al. (1997) adopt the mass of the 
hydrogen-free core at the start of the first thermal pulse 
as the mass of the resulting white dwarf. 

The implication of this dichotomy in the mass of the 
progenitor star for the white dwarf carbon-oxygen com- 
position is illustrated in Fig. [3] We show the inner car- 
bon and oxygen abundance distribution for the 0.63 Mq 
white dwarf resulting from the two different progenitors 
discussed previously. For illustrative purposes, chemi- 
cal rehomogenization by Rayleigh- Taylor instabilities has 
not been considered in this particular example. The thick 
lines display the chemical profile for the progenitor star 
with 3 Mq characterized before the first thermal pulse by 
a hydrogen-free core of 0.63 Mq. The thin lines show the 
chemical profile for the 2.25 Mq progenitor which leads to 
the same white dwarf mass but after evolution has pro- 
ceeded through the thermally-pulsing AGB phase (see 
Fig (21). Note the different chemical profiles expected in 
both cases. In particular, the central oxygen abundance 
may be underestimated by about 15% should the white 
dwarf mass is assumed to be the hydrogen-free core mass 
right before the first thermal pulse. Note that, however, 
this variation is an upper limit, and it would be less for 
other white-dwarf masses. Clearly, the initial-final mass 
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Fig. 5. — Abundance distribution of hydrogen, helium and car- 
bon in terms of the outer mass fraction for the 0.878, 0.609 and 
0.525 Mq white dwarf models at two selected stages just after the 
maximum effective temperature point and near the beginning of 
the ZZ Ceti regime (thin and thick lines, respectively). 

relationship is an aspect that has to be considered in the 
problem of the carbon-oxygen composition expected in 
a white dwarf, as well as in attempts at constraining, 
from pulsational inferences of variable white dwarfs, the 
mixing processes and the efficiency of the ^^C(a,7)^®0 
reaction rate in the core of helium burning stars. Also, it 
is important to realize that a larger oxygen abundance in- 
creases the cooling rate of the white dwarf because of the 
lower heat capacity and because an initial larger oxygen 
abundance reduces the energy release by phase separa- 
tion on cyrstallization (Isern et al. 2000). 

4. THE INTERNAL CHEMICAL PROFILES 

The carbon-oxygen stratification for some selected 
models is shown in Figs. |3] and |31 The shape of the 
chemical profiles before rehomogenization is given in Fig. 
[31 Easily recognizable are the flat chemical profiles in 
the inner part of the core left by convection during the 
core helium burning, the bumps resulting from the inclu- 
sion of extra-mixing episodes beyond the fully convective 
core, and the signatures of the outward-moving helium- 
burning shell after the end of core helium burning. Be- 
cause of the larger temperatures in the helium burning 
shell, the oxygen abundance decreases in the outer re- 
gions of the carbon-oxygen core. 

The expected chemical profiles of some of our white 
dwarf models when evolution has proceeded to the do- 
main of the ZZ Ceti stars are displayed in Fig. 2] At this 
stage, chemical rehomogenization by Rayleigh- Taylor in- 
stabilities has already smeared out the bumps in the in- 
ner profiles, leading to quite extended flat chemical pro- 
files. Note the dependence of both the core chemical 



abundances and the location of the chemical transitions 
on the stellar mass. Pulsation periods in white dwarf 
models are very sensitive to the shapes and locations of 
the chemical transions zones (see section 5). This empha- 
sizes the need for a detailed knowledge of the progenitor 
history for a realistic treatment of white dwarf evolution 
and pulsations. The chemical profiles in the outermost 
regions, resulting from prior mixing and burning events 
during the thermally-pulsing AGB phase, are markedly 
modified by the diffusion processes acting during white 
dwarf evolution, particularly in the case of more mas- 
sive models, where chemical diffusion and gravitational 
settling are notably more efficient. This can be better 
appreciated in Fig. [S] where the chemical abundance 
distribution for white dwarf models of different stellar 
masses is depicted in terms of the outer mass fraction. 
These plots put special emphasis in the outer regions of 
the model. In this figure and for each stellar mass, thin 
lines show the chemical abundance distribution at early 
stages of white dwarf evolution when diffusion has not 
had time to act. The signature of the evolution through 
the thermally-pulsing AGB stage on the chemical pro- 
file, particularly the formation of the helium-rich buffer 
and the underlying intershell region rich in helium and 
carbon — built up during the mixing episode at the last 
AGB thermal pulse — are easily visible. The presence 
of carbon in the intershell region stems from the short- 
lived convective mixing that has driven the carbon-rich 
zone upward during the peak of the last helium pulse on 
the AGB. Thick lines depict the situation at advanced 
stages, near the ZZ Ceti instability strip, when element 
diffusion has strongly modified the chemical abundance 
distribution and resulted in the formation of a thick pure 
hydrogen envelope plus an extended inner tail. In the 
more massive models, chemical diffusion leads to a sig- 
nificant amount of carbon in the helium buffer zone. Also 
near-discontinuities in the initial abundance distribution 
are smoothed out considerably by diffusion. 

An important fact to note in Fig. [5] is the dependence 
on the stellar mass of the outer layer chemical stratifica- 
tion expected in ZZ Ceti stars. Indeed, for the more mas- 
sive models, diffusion has strongly modified the chemical 
abundance distribution, eroding the intershell region by 
the time evolution has reached the domain of the ZZ 
Ceti instability strip. This is in contrast with the sit- 
uation encountered in our less massive models, where 
the intershell region is not removed by diffusion. This is 
because element diffusion is less efficient in less massive 
models and also because the intershell is thicker, with the 
subsequent longer diffusion timescales. Regarding white 
dwarf asteroseismology, these are not minor issues, since 
the presence of a double-layered structure in the helium- 
rich layers is expected to affect the theoretical 5— mode 
period spectrum of ZZ Ceti stars, as it does for pulsating 
DB (He- rich) white dwarf stars (Metcalfe et al. 2003). 
However, we mention that the thickness of the intershell 
region also depends somewhat on the number of thermal 
pulses during the AGB experienced by the progenitor 
star. 

Finally, we have explored whether the shape of the 
He/H chemical interface has a dependence on the thick- 
ness of the H envelope. This is an important issue, be- 
cause predictions of the exact value of the H envelope 
mass are tied to the precise mass loss history along the 



previous AGB and post-AGB phase, and particularly to 
the occurrence of late thermal pulses. We have gen- 
erated additional, artificial sequences with H envelopes 
much thinner than those of our models. These sequences 
were created at high luminosities from a hot model with 
M* = O.6O9M0. This artificial procedure took place at 
luminosities high enough as to ensure that the models be- 
come physically sound at stages far before the domain of 
the ZZ Ceti stars. We found that diffusion rapidly leads 
to pure hydrogen envelopes, but the shape of the He/H 
chemical interfaces by the time evolution has proceeded 
to the ZZ Ceti stage, is almost the same independently 
of the thickness of the H envelope. 

5. PULSATION PROPERTIES: COMPARISON WITH 
PREVIOUS CALCULATIONS 

White dwarf asteroseismology is sensitive to the pre- 
cise shape of the internal chemical profiles. The entire 
g—mode period spectrum and mode-trapping properties 
of pulsating white dwarfs are very sensitive to the fine 
details of the chemical profiles of both the core and the 
envelope of the star. This extreme sensitivity has been 
exploited with some success in several pulsation stud- 
ies to infer the core chemical structure — e.g., Metcalfe 
(2003) for the case of pulsating DB white dwarfs — and 
the thickness of the He and H envelopes — e.g., Bradley 
(1998, 2001) for DAV stars. 

In this section we perform a comparison between the 
pulsation properties derived from our new chemical pro- 
files and those based on the most widely used chemi- 
cal profiles in existing asteroseismological studies. To 
assess the adiabatic pulsation properties of our white 
dwarf models we employ the numerical code described 
in Corsico & Althaus (2006). We refer the reader to that 
paper for details. With the aim of simplifying our anal- 
ysis, we elect a template DA white dwarf model with 
JVh = 0.6096 Mq, Toff ~ 12 000 K, and a thick hydrogen 
envelope (Mh ~ lO^'^Af,). This is a canonical model of 
a DAV star with an average mass located in the mid- 
dle of the observed ZZ Ceti instability strip. In Fig. 
[niwe depict the internal chemical profile (upper panel), 
the Ledoux term (middle panel), and the logarithm of 
the squared Brunt- Vaisala and Lamb frequencies (lower 
panel) for dipole modes corresponding to our template 
model. The Ledoux term and the Brunt- Vaisala fre- 
quency are computed as in Corsico & Althaus (2006). 
Our model is characterized by three chemical transition 
regions: a double chemical interface of oxygen and car- 
bon located at the core region (0.2 < — logg < 0.8), a 
triple chemical interface of oxygen, carbon, and helium 
located at 1.2 < — logq < 2.3, and finally, a double chem- 
ical interface of He/H located at 3.0 < — logq < 4.5. The 
core chemical profile is typical of situations in which ex- 
tra mixing episodes beyond the fully convective core (like 
overshooting) during the core He burning phase are al- 
lowed to operate. The smoothness of the He/H chemical 
interface, on the other hand, is the result of the time- 
dependent element diffusion processes. The existence of 
these three chemical interfaces induces the "bumps" in 
the profile of the Brunt- Vaisala frequency. The num- 
ber of these bumps, as well as their heights and widths, 
strongly affect the whole structure of the pulsation spec- 
trum of the star. 

In Fig. [7] we show the situation in which our tem- 
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Fig. 6. — The internal chemical profile (upper panel), the Ledoux 
term B (middle panel), and the logarithm of the squared Brunt- 
Vaisala (A'') and Lamb (Li) frequencies (lower panel) in terms of 
the outer mass fraction (— log q, where q = 1 — Mr/Mt) for dipole 
modes (^ = 1) corresponding to a DA white dwarf model with 
A/* = 0.6096 Mq and T<,ff ~ 12 000 K. 

plate model is characterized by a ramp-like core chem- 
ical structure of the kind used by Bradley (1996, 1998, 
2001) and more recently by Bischoff-Kim (2009). In this 
case the core chemical profile is not the result of stellar 
evolution calculations, but parameterized. This kind of 
chemical profiles has been widely employed in asteroseis- 
mology of white dwarfs because it allows a full explo- 
ration of the parameter space regarding the shape of the 
chemical abundance profiles in the core. The parame- 
ters are the central oxygen abundance (Xq), the coordi- 
nate at which Xq starts to drop, and the coordinate at 
which Xo drops to zero — see, e.g., Bischoff-Kim et al. 
(2008). The overall shape of the core chemical profiles 
displayed in Fig. [7] is substantially simpler than that of 
the chemical profiles characterizing the model depicted 
in Fig. [6l Note the presence of a chemical transition re- 
gion of C and He. This is at variance with the chemical 
profiles produced by LPCDDE, which are characterized by 
a triple transition region of oxygen, carbon, and helium 
(see Fig. [S]). The shape of the C/He chemical interface 
in this model is set by diffusion parameters, which were 
chosen to match the pulsation periods of G 117— B15A 
(Bischoff-Kim et al. 2008). It is worth noting that this 
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Fig. 7. — Same as in Fig 
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|6] but for the case of ramp-like core 
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Fig. 8. — Same as in Fig. |6] but for the case of Salaris-like core 
chemical profiles. 



transition region produces the most prominent bump in 
tlie Brunt- Vaisala frequency (see the Ledoux term B in 
the central panel of Fig. [T]). We mention that the pres- 
ence of a thick pure carbon buffer like that assumed in 
these models is not expected from stellar evolution cal- 
culations. Finally, the He/H chemical transition region 
has been obtained by assuming equilibrium diffusion (Ar- 
coragi & Fontaine 1980), but without the trace element 
approximation — see Bischoff-Kim et al. (2008) for de- 
tails. The shape of this interface is very similar to one 
based on time-dependent diffusion. It is worth noting, 
however, that the slight differences in the thickness and 
steepness of this chemical interface between the LPCODE 
model and the ramp-like model lead to a non-negligible 
contribution to the differences found in the period spac- 
ing (and the period themselves) of the models (see later 
in this section). 

In Fig. [5] we display the case of a template model 
characterized by a Salaris-like core chemical structure. 
This kind of core chemical profiles was employed first 
by Corsico et al. (2001, 2002a) and Benvenuto et 
al. (2002a, b), and more recently in Bischoff-Kim et al. 
(2008). Actually, the core chemical profiles displayed in 
Fig. m and those used in Bischoff-Kim et al. (2008) 
are a close adaptation of the original chemical profiles of 
Salaris et al. (1997). Except for the core region, the rest 



of the chemical profiles in this model are the same as in 
the model depicted in Fig. [71 At variance with the tem- 
plate models described before, in this case there are two 
core chemical interfaces of oxygen and carbon instead of 
just one. This leads to four bumps in the Brunt- Vaisala 
frequency, as can be seen in the lower panel in Fig. [S] 

Figures IH [3 and [8] clearly reveal the profound differ- 
ences existing between the chemical profiles and the run 
of the Brunt- Vaisala frequency of the model generated 
with LPCODE, i.e., by considering the full evolution of 
progenitor stars, and the two template models that in- 
corporate the most widely used chemical profiles in past 
and current asteroseismological studies of ZZ Ceti stars. 
The differences are particularly noteworthy in the core 
chemical structure. 

In what follows, we compare the pulsation properties of 
our template model with the models having core ramp- 
like and Salaris-like chemical profiles. The £ — 1 asymp- 
totic period spacing is largest for the LPCODE model 
(45.38 s), followed by the Salaris-like core chemical pro- 
file model (44.17 s) and by the ramp-like model (43.32 
s). Since these models have the same stellar mass and 
effective temperature, the period and asymptotic period 
spacing differences are exclusively due to the differences 
in the chemical profiles at the core and the envelope of 
the three models. In particular, the subtle differences ex- 
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Fig. 9. — Comparison between the template model generated 
with LPCODE (Fig. [G} and the template model with ramp-like core 
chemical profiles (Fig. ^. The upper panel shows the differences 
in the k identification, and the lower panel depicts the differences 
between the matched periods. 

isting in the shape of the He/H interface resuhing from 
time-dependent diffusion and that obtained from equihb- 
rium diffusion give rise to a non-negligible contribution 
to the difference in the asymptotic period spacing. 

If we now turn the argument around and imagine doing 
asteroseismological fits where we fix the chemical profiles 
to those found using LPCODE and allow the mass and ef- 
fective temperature of the models to vary to give us a 
pulsation spectrum that matches that of an observed pul- 
sating white dwarf, the asymptotic period spacing of the 
models will influence the mass and effective temperature 
of the best fits models. As it is well known, the hotter 
and more massive models have smaller asymptotic period 
spacings (for modes that are not strongly trapped) . Since 
the asymptotic period spacing is larger for the LPCODE 
models, the best fit models would have to have larger 
mass and effective temperature to match a given ob- 
served asymptotic period spacing. As a result, we would 
expect to find asteroseismological fits that are more mas- 
sive and hotter than current fits. This effect should 
mainly be observed for rich white dwarf pulsators, where 
we have a wealth of higher k (asymptotically spaced) 
modes to fit. For G 117— B15A, for instance, we cannot 
draw any conclusions from the asymptotic period spac- 
ings alone, as the only 3 modes observed have low radial 
overtone and are strongly trapped (Bischoff-Kim et al. 
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Fig. 10. — Same as in Fig. |9l but for the comparison between 
the template model generated with LPCODE (Fig. |6ll and the 
template model with Salaris-like core chemical profiles (Fig. \8^. 



2008). 

The higher period spacing for the LPCODE model leads 
to a drift to longer and longer periods as we work down 
the list of periods toward higher k modes. For instance, 
the cumulated effect of the 2.06 s difference between the 
LPCODE and the ramp-like models results in higher k pe- 
riods to differ by as much as 65 seconds (for A; '^ 30 if 
£ = 1 and fc ~ 55 if i' = 2). The practical result in as- 
teroseismological studies of such a drift to higher periods 
would be to lead to different k identifications of modes. 
In Figs. HI and [10] we show how asteroseismological fits 
may be affected. In the upper panels, we show the differ- 
ences in k identifications and in the bottom panels, how 
the period of the matched modes differ. Even allowing 
the k identifications to change to find the best match be- 
tween the periods of the two models compared, we still 
find that some individual periods differ by as much as 
^ 25 seconds. 

We conclude our analysis by examining the forward 
period spacing (Allfc = Ilfc+i — 11^) and the kinetic en- 
ergy (£'kin) of the three template models (Figure [TT\ . 
The kinetic energy is computed according to Corsico & 
Althaus (2006). The horizontal lines in red correspond 
to the asymptotic period spacing. As it is well known, 
the presence of density gradients in the chemical transi- 
tion regions in the interior of a white dwarf causes mode 
trapping (Winget et al. 1981; Brassard et al. 1992a; 
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of ^ = 1 modes in terms of the pulsation periods. The horizontal lines in the upper panels display the asymptotic period spacing. 



Corsico et al. 2002b). The signature of mode trapping on 
the period-spacing distribution is the presence of strong 
minima in a Allfc — 11^ diagram, in contrast to the sit- 
uation in which the star is chemically homogeneous — 
see, for instance, Corsico & Benvenuto (2002). Due to 
the presence of several chemical interfaces in the tem- 
plate models, we expect to find clear signatures of mode 
trapping. 

For the model generated with LPCODE, the Allfc dis- 
tribution exhibits two primary minima and several sec- 
ondary minima, as can be seen in the upper left panel of 
Figure [TT] In the case of the model with ramp-like core 
chemical profiles (upper central panel) there are also pri- 
mary and secondary minima, but the contrast amongst 
them is much less pronounced than in the case of the 
LPCDDE model. Finally, in the case of the model with 
Salaris-like core chemical profiles (upper right panel) 
there is no clear distinction between primary and sec- 
ondary minima of Allfc. That is, all the minima are very 
similar. 

In spite of the complexity of the mode-trapping pat- 
tern exhibited by the template models, it is possible to 
make some rough inferences by examining the values of 
the kinetic energy of the modes (lower panels of Fig. 
[TT|) . A close inspection of the plots reveals that each 
minimum of Allfe is associated with a maximum in -Ekin 
of a mode with the same radial order k or differing in 
Afc = ±1. Modes corresponding to local maxima in 
ii^kin are modes partially confined to the core regions be- 
low the 0/C and/or the C/He chemical interfaces (the 
0/C/He transition in the case of the LPCDDE model), that 
is, modes with amplitudes relatively large even in very 
deep layers of the model. These modes are potentially 
useful from an asteroseismological point of view. The 
remaining modes (which have neither maxima nor min- 
ima kinetic energy values) are much less sensitive to the 



presence of the chemical transition regions, and so, they 
are of minor relevance for asteroseismolgy 

We conclude that the pulsation properties of DA white 
dwarf models that incorporate our new chemical profiles 
for the core and envelope substantially differ from those 
of models having the most commonly used chemical pro- 
files. The important issue to be addressed now is to asses 
the impact of our new chemical profiles on asteroseismo- 
logical period-to-period fits of real DAV stars. This step 
is beyond the scope of the present work, and we defer it 
to future papers. 

6. SUMMARY AND CONCLUSIONS 

In this paper we computed new chemical profiles for 
the core and envelope of white dwarfs appropriate for 
pulsational studies of ZZ Ceti stars. These profiles were 
derived from the full and complete evolution of progen- 
itor stars from the zero age main sequence, through the 
thermally-pulsing and mass-loss phases on the asymp- 
totic giant branch (AGB). These new profiles are in- 
tcnted for asteroseismological studies of ZZ Ceti stars 
that require realistic chemical profiles throughout the 
white dwarf interiors. In deriving the new chemical pro- 
files, we employed the LPCDDE evolutionary code, based 
on detailed and updated constitutive physics. Extra- 
mixing episodes during central hydrogen and helium 
burning, time-dependent element diffusion during the 
white dwarf stage and chemical rehomogenization of the 
inner carbon-oxygen composition by Ray leigh- Taylor in- 
stabilities were considered. The metallicity of progenitor 
stars is Z ~ 0.01. 

We discussed at some length the importance of the 
initial-final mass relationship for the white dwarf carbon- 
oxygen composition. A reduction of the efficiency of 
extra-mixing episodes during the thermally-pulsing AGB 
phase, supported by different pieces of theoretical and 
observational evidence, yields a gradual increase of the 
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hydrogen-free core mass as evolution proceeds during this 
phase. As a result, the initial-final mass relationship by 
the end of the thermally-pulsing AGB is markedly differ- 
ent from that resulting from considering the mass of the 
hydrogen free core right before the first thermal pulse. 
We found that this issue has implications for the carbon- 
oxygen composition expected in a white dwarf. In par- 
ticular, the central oxygen abundance may be underesti- 
mated by about 15% if we assume the white dwarf mass 
to be the hydrogen-free core mass before the first thermal 
pulse. We also discuss the importance of the computa- 
tion of the thermally-pulsing AGB and element diffusion 
for the chemical profiles expected in the outermost lay- 
ers of ZZ Ceti stars. In this sense, we found a strong 
dependence of the outer layer chemical stratification on 
the stellar mass. In less massive models, the intershell 
region rich in helium and carbon — which is built during 
the mixing episode at the last AGB thermal pulse — is 
not removed by diffusion by the time the ZZ Ceti stage 
is reached. 

Finally, we performed adiabatic pulsation computa- 
tions and discussed the implications of our new chem- 
ical profiles for the pulsational properties of ZZ Ceti 



stars. We found that the whole g— mode period spec- 
trum and the mode-trapping properties of these pulsat- 
ing white dwarfs as derived from our new chemical pro- 
files are substantially different from those based on the 
most widely used chemical profiles in existing asteroseis- 
mological studies. 

We expect the best fit parameters of asteroseismolog- 
ical studies using the LPCODE chemical profiles to differ 
significantly from those found in studies made so far. 
Further studies will show in what way. Will we solve 
the high mass problem with G117— B15A and Salaris- 
like core profiles (Bischoff-Kim et al. 2008) or find thicker 
hydrogen layers in asteroseismological fits, more in line 
with stellar evolution calculations (Castanheira & Kepler 
2008)? 
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